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Abstract 

We study the process of thermal activation mediated by sphaleron transitions by 
analyzing the real-time dynamics of the decay out of equilibrium in a 1 + 1 dimensional 
field theory with a metastable state. The situation considered is that of a rapid super- 
cooling in which the system is trapped in a metastable state at a temperature larger 
than the mass of the quanta, but smaller than the energy to create a critical droplet. 
The initial density matrix is evolved in time and the nucleation rate (probability cur- 
rent at the saddle point) is computed. The nucleation rate is time dependent, vanishing 
at early times, reaching a maximum at a time t ~ 1/m with m the mass of quanta in 
the metastable state, and decreasing at long times as a consequence of unitarity. An 
estimate for the average number of particles of "true vacuum" produced as a function 
of time during the nucleation process is obtained. 
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1 Introduction and Motivation 



Thermal activation plays a very important role in the dynamics of evolution out of equilib- 
rium of thermodynamically metastable states. Such a situation arises for example in first 
order phase transitions which may have ocurred in the early universe|l|, 0, ^, where false 
vacuum decay provides a mechanism for ending inflation. More recently thermal activa- 
tion has been conjectured to be the mechanism responsible for unsuppressed baryon number 
violation at high temperatures |^. Within the context of baryon number violation, 
Kuzmin, Rubakov and Shaposnikov0, argued that in non-Abelian gauge theories topolog- 
ically different vacua (characterized by different Chern-Simons number) are separated by 
energy barriers and that there are configurations (in functional space) that are responsible 
for over the barrier transitions from one vacuum to another. Klinkhamer and Manton^ 
suggested that these over-the-barrier transitions correspond to static configurations that ex- 
tremize the energy functional with at least one unstable mode. These authors named this 
configuration the "sphaleron" and provided a topological proof for the existence of such con- 
figurations in the Standard Model. In non-Abelian gauge theories this configuration carries 
half a unit of winding number, and corresponds to a saddle point in the space of functions. 
These transitions are necessarily accompanied by a change in baryon number [§, ^ p!T[] . 



Within the context of first order phase transitions at finite temperature, LindefT^ argued 
that the relevant configuration responsible for the decay of the metastable state is an 0(3) 
(static) "bounce", which is in a sense the high temperature limit of the zero temperature 
bounce that describes a tunneling event |ll3| between the metastable and stable phases. 

An important quantity is the transition rate (per unit volume), as this quantity deter- 
mines the relevant time scales for the completion of first order phase transitions in infla- 
tionary cosmologies and for the rate of baryon number violation in gauge theories at finite 
temperature. 

The usual calculation of the rate follows hanger's original theory of homogeneous nucle- 



ation in classical statistical mechanics |Tj], and Affleck's|T^ generalization to the quantum 
case. 

hanger's treatment is based on the computation of the probability current at the saddle 
point corresponding to a steady state solution of an appropriate Fokker-Planck equation. It 
is assumed in this treatment, and clearly spelled out in hanger's calculation, that the steady 
state condition corresponds to the situation in which the metastable phase is replenished 
at the same rate at which probability is flowing out across the saddle. There is a "source" 
of probability on the metastable well and a corresponding "sink" on the stable side. The 
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saddle point configuration corresponds to a "critical droplet" . The probability flows from the 
metastable to the stable phase along the unstable direction in functional space. Under the 
steady-state assumption the rate is calculated in equilibrium and identified to be proportional 
to the imaginary part of an analytically continued free energy. The same assumptions are 
implicit in Affleck's approach in which the transition rate is computed as an equilibrium 
average |ll5[|. As a result of the steady-state assumption, the rate is independent of time and 
basically determined by an Arrhenius (thermal activation) equilibrium Boltzmann factor [|16| . 



hanger's treatment of nucleation has been recently generalized to relativistic fleld theory 
(under basically the same assumption of a steady state) by Csernai and Kapusta^]. 

The "critical droplets" in hanger's theory of homogeneous nucleation are thus identifled 
with the "sphalerons" in the broader sense studied by Manton and Samols||T8[], i.e, static 
fleld conflgurations with (at least) one unstable mode. These are saddle points of the energy 
functional corresponding to fleld conflgurations "sitting" at the top of the energy barrier (in 
functional space). Numerical simulations of sphaleron transitions have been carried out by 
Hellmund and Kripfganzfl^ and Ambjorn et. al.[^. These authors integrate numerically 
the classical equations of motion by starting from a sphaleron conflguration. Recently a 
numerical study of the process of thermal activation assuming a simple Langevin description 
has been reported . 



The motivation of this article is to provide a different approach to the study of the decay 
of a metastable state by considering the process of thermal activation through "sphaleron 
transitions" as a non equibrium evolution in real time. 

We consider the situation of a 1 + 1 dimensional scalar fleld theory in which a metastable 
state has been formed after a period of rapid supercooling and decays via the formation of 
"critical droplets" and their subsequent evolution. We propose to describe the supercooled 
metastable state in terms of a functional density matrix corresponding to a thermal distri- 
bution of free-fleld quanta (harmonic oscillator states) of the false vacuum. The temperature 
of the ensemble is assumed to be much larger than the mass of the quanta (m) but smaller 
than the energy to create a sphaleron conflguration. The "sphaleron" determines an unsta- 
ble direction in functional space along which probability flows. We solve for the real time 
evolution of the initial density matrix along this unstable direction. The transition rate is 
identifled with the probability current flowing across the saddle from the metastable to the 
stable state. This approach provides a description of the decay process out of equilibrium 
from an initial state. We flnd that the transition rate is time dependent and that it strongly 
depends on the initial state. 

In section II we study the properties of the "sphaleron" and quantize the theory around 
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this static semiclassical configuration. Section III presents a detailed analysis of the time 
evolution of the field theoretical density matrix and the probability current along the unsta- 
ble direction on the saddle point is evaluated. Following Langer's original treatment, this 
probability current is identified with the transition rate. In this section we also mention 
some recent experiments on nucleation in classical fluids that report a time dependent rate, 
and dependence on initial conditions. 

In section IV we present an approximate calculation for the production of particles of the 
true vacuum during the decay process. Section V summarizes our conclusions, establishes 
the range of validity of our results, and explores potential implications. 

Two appendices contain relevant technical details. 



2 The Sphaleron and Thermal Activation: 

We consider a scalar field theory in 1 + 1 space-time dimensions with Hamiltonian 



H 



'd^{x) 



dx 



+ Vmx)) 



(2) 



The potential is depicted in figure 1(a). The point $* is parametrized by the mass m of 
harmonic oscillations around the minimum $_ as 



1 



m 



(3) 



In two space time dimensions, the scalar field is dimensionless and the coupling constant has 
dimensions of {massy. 

The "sphaleron" is defined as a static field configuration that corresponds to an extremum 
of the energy functional with at least one unstable direction in functional space. In our case, 
a stationary solution of the equations of motion (extremum of the energy functional) satisfies 

_ d'^x) dVi<!>) ^ 

dx^ 9$ ^ ' 

This equation resembles that of a particle moving in "time" (labeled by x) down a potential 
—V. The solution that starts at $_ at x ^ — oo and returns to as a; — oo is easily found 
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where Xo is an arbitrary integration constant and reflects the translational invariance of the 
equations of motion. This solution corresponds to a kink-antikink pair or "droplet" with 
a "radius" 2so/m, and is similar to a "polaron solution" found in quasi-one dimensional 
polymers ||23|]. This solution and its stability was also studied by Baltar et. al.p4| 

The dimensionless ratio e is a measure of the depth of the global minimum ($+). As 
e —* 1 the minima of the potential become degenerate and So — oo. The situation e ~ 
1; So ^ 1 corresponds to a "thin- wall" droplet for which the radius is much larger than the 
wall thickness or "skin" of the one-dimensional droplet ^ = 2/m, as depicted in figure 2. For 
a "thin-wall" droplet we find the energy of the sphaleron configuration to be 

.3 
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where we have shown that the corrections are exponentially small in the "thin- wall" limit. 

It is convenient to expand the field and its canonical momentum around the sphaleron 
configuration as 



<l>(x) 
U(x) 



Xq) 



-ih6i 



(9) 
(10) 

(11) 



with the fi{x — Xo) being the eigenf unctions of the operator of quadratic fluctuations around 
the sphaleron configuration 



dx' 



+ V"{<^sph{x-Xo)) 



fl{x - Xo) = UJi fi{x - Xo) 



(12) 



chosen to be real and orthonormal (in a volume L). We have explicitly written the dependence 
of the eigenfunctions on Xo as a consequence of translational invariance, and indicated that 
the operators (pi, tti depend parametrically on Xo- This will become important later when we 
introduce collective coordinates to treat translational invariance. 

Translational invariance guarantees that there is a zero frequency mode whose normalized 
eigenf unction is given by 

1 d^sph{x - Xo) 



foix - Xo 



sph 



dx 



(13) 
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This eigenfunction is parity odd and has one node, and in the thin-wall limit it is identified 
as the antisymmetric linear combination of the zero modes for kink and anti-kink. There is 
an orthogonal symmetric combination that is parity even and nodeless corresponding to the 
orthogonal linear combination, which is the eigenfunction with the lowest eigenvalue. In the 
thin-wall limit we find it to be 

J. _ m d^sphjx - Xo\ gp) 

This mode has a negative eigenvalue ijJ^_^ = —Q^ [Q > 0), in the "thin- wall" limit this 
eigenvalue is exponentially small, determined by the overlap integral of the kink-antikink 
zero modes 0{e~'^^°). In this approximation, and based on the known results on 

bound-states of the one kink case, we conclude that the rest of the spectrum is positive 
definite. 

The Hamiltonian in this basis, to be referred to as the "sphaleron basis" , becomes 



H = Esph + H^ + Hi (15) 
= f + 4^-^0-1 + ^EK + a;f0.] (16) 

where Hj, has terms cubic and quartic in terms of the 0; and we have suppressed the upper 
index [x^ in both the tt/ ; (f)i to avoid cluttering of notation, but keeping in mind that these 
operators depend parametrically on {xo). 

The nature of the instability represented by the mode Z = — 1 is physically clear. For 
small amplitudes of 0o ; 0-i the field may be written as 

^{x) ^ ^sph[x-{xo + 5xo)\{so + 5so)]+Y,h{^-^o)(t>i (17) 

i>i 

5xo = _^2= ; 5so = — ^=0_i (18) 

Whereas 0o represents a translation of the center of mass of the sphaleron configuration, 0_i 
represents an expansion or contraction of the radius of the droplet. Thus we identify the 
collective coordinate that describes the radius of the droplet (in this approximation) as 

Tfi 

s = So+^^=(p-i (19) 
Egph 

The coordinate measures the departure from the critical droplet (for small amplitudes 
of this coordinate). 
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When the droplet expands, the configuration gains volume energy because it is sampling 
a larger region in space with a lower energy density. The gain in volume energy grows linearly 
(asymptotically) in one space dimension. 

The gradient terms, giving rise to a surface energy of the droplet, saturate at about twice 
the kink mass in one space dimension. To see this clearly, let us define a field configuration 
(droplet) parametrized by the radius s, as 



and use s as a parameter. The (classical) energy density as a function of this parameter 
£{x, s) is depicted in figure 3 (a, h, c) for s < So ; s = So ; s > So respectively, for the case of 
the potential with a metastable minimum (figure 1(a)). 

The total (classical) energy E[s] is depicted in figure (4) for the potentials shown in figure 
(1). The maximum of E[s] for the case of a potential with a metastable minimum (figure 5) 
is given by So given by equation (|^), as this is the value that corresponds to the solution of 
the equation of motion for the static configuration, that is, the "sphaleron". The value So 
corresponds to a "critical droplet", for s < So the droplet will shrink, as the cost in elastic 
surface term is greater than the gain in volume energy. A droplet with s = Sq is in unstable 
equilibrium, whereas for s > So the droplet will grow as the gain in volume energy offsets 
the cost of elastic surface (wall) energy. 

The sphaleron configuration corresponds to a saddle point in functional space; there is 
one unstable direction corresponding to the dilation or contraction of the droplet, one fiat 
direction corresponding to translational invariance, and the remaining (infinite) directions 
are all stable, with the energy increasing quadratically for small amplitudes away from the 
saddle. 

This unstable direction in functional space plays a fundamental role: consider an initial 
quantum mechanical state localized in the metastable minimum, for example either a ground 
state wave functional or a thermal density matrix for the quadratic potential centered on 
the metastable minimum (see figure 1(a)). If the probability for finding a droplet with a size 
greater than the critical size in the initial state is non-zero, the time evolution of this state will 
inexorably tend to spread the state in such a way that field configurations which will sample 
the global minimum in larger regions in space will acquire larger probabilities as time evolves. 
This translates into the notion that the time evolved state will give rise to a probability 
current at the saddle along the unstable direction towards field configurations corresponding 
to growing droplets. The "decay rate" of the initial state (or ensemble) will be determined 
by the total probability current along this direction in functional space. Furthermore, by 




(20) 
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translational invariance, these droplets appear with equal probability at all points in space, 
and the total current passing through the saddle point along the unstable direction will be 
proportional to the volume. 

At finite temperature, there will be a non-zero probability of finding a critical (and larger) 
droplet in the initial ensemble and the instability will take the initial ensemble away from 
equilibrium. 

The "decay" of the initial ensemble will thus correspond to a process of thermal activation, 
in which the initial state is driven "over the barrier" in functional space. This barrier should 
not be confused with the hump in the scalar potential, but it must be identified with the 
maximum of E[s] in figure 5 for the metastable case, that is to say, the energy of a critical 
droplet. 

We now study this process of thermal activation in real time by following the time evo- 
lution of an initially prepared density matrix. 

3 Real Time evolution of Initial Ensembles: 

We consider the situation in which there is a very rapid supercooling from a high temperature 
phase in thermal equilibrium to a situation in which the system gets trapped in a metastable 
state. This corresponds to the case in which the initial ensemble is described by a density 
matrix in thermal equilibrium, but "centered" at the metastable minimum of the potential. 
We approximate this initial density matrix as that of a free field theory centered at the 
metastable minimum and with particles with mass m determined by the second derivative 
of the potential at In the Schroedinger representation the density matrix elements are 
given by 




A^exp [-2J J [(*(^) - - y)i.^i.y) - ^-)+ 

-2($(a;) - ^.)K2{x - y){^\y) - $_)]} 




(21) 




y) 




(22) 




y) 



N 




(23) 



(24) 
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We have chosen the normahzation factor in such a way that 

Trp = j V^p{^,^) = 1 (25) 

With this normahzation, the quantity p($, $) has the interpretation of a probabihty density 
in functional space. 

In terms of the Fourier components of the field $, the above density matrix is easily iden- 
tified with a density matrix in the coordinate representation for a collection of independent 
harmonic oscillators with frequencies uj^- 

In particular, a meaningful question is: what is the probability (density) to find the 
sphaleron configuration in the initial ensemble?. This quantity is easily calculated in the 
high temperature limit where the Fourier transform of ($sp/i(x) — $_) is non-negligible only 
for momenta k ^ m. Thus in the convolution of the sphaleron configuration and the kernels, 
for temperatures T ^ m, the kernels in the density matrix may be approximated by 



Ki{k)-K2{k) = ^tanh 



2kBT 



Iksi 



after some straightforward algebra we find in the "thin-wall" limit 



p(<l>spft,$spft) iVexp 



2kBT 



1 + 12sJ 



(27) 



with So the radius of the droplet. 

Then we see that the probability of finding the "sphaleron" configuration in the initial 
ensemble is very different fTom exp[—Esph/kBT]. The latter corresponds to a density matrix 
of the form po = ^w[—{Esph + Hq)/kBT], with Hg the Hamiltonian (|TBp. This density matrix 
represents an ensemble "centered" at the sphaleron configuration. 

We can also compute the "average radius" of the droplets in the initial state. We find 
(see appendix B) 

^ --\/^ (28) 
(.) ^ so-l (29) 

in the thin wall limit Sq ^ 1 most of the "droplets" in the initial ensemble are slightly 
smaller than the critical droplet. This is clearly in agreement with the condition that our 
initial ensemble is describing the metastable phase and the field configurations that sample 
large regions of the stable phase are very rare. 
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The time evolution of the initial density matrix (^) is easily obtained in the quadratic 
approximation for ([T6|) , that is setting Hi = 0. The reason is that both the initial density 
matrix and the Hamiltonian (around the sphaleron configuration) are quadratic, thus at all 
times the density matrix will be of the gaussian form. By expanding the field in the sphaleron 
basis and defining 

$_-$,p/,(x-Xo) = ^/i(x-Xo)0z ; 00 = (30) 

I 

the initial density matrix becomes 

p(0,0',t = O) = iVexp|-i[(iri),,;,(r/,r/;, + r/Wr)-2(ir2)M'Wr]} (31) 
(i^i)M' = Jd^JdyJ ^e-''^^-y^K,ik)fi{x - Xo)fv{y - Xo) (32) 

(i^2)M' = Jd^JdyJ ^e-''^'-y^K2{k)Mx - Xo)My - xo) (33) 
VI = (0/ - 0z) ; V'l = (0'^ - 00 (34) 

Thus we see that in the sphaleron basis, the different modes are mixed in the initial density 
matrix. 

The time evolution of the density matrix is formally determined by the Liouville equation 

^f^^p(Ml2ll = [H,pi<!>,<^',t)] (35) 

whose formal solution is 

p($, t) = e-^"'p{<^, 0)e^^* (36) 

In the Schroedinger representation with Il{x) = —ih6/6^{x) it is straightforward to see 
that the functional probability density obeys a continuity equation 



dt J 5^{x 

In the "sphaleron basis" , the above continuity equation becomes 

JAM = ^(^r^)p^M\t)U-_, (40) 
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we have suppressed the label (xo) , but we must keep in mind (and this will become relevant 
later) that the (pi depend on it (though the (pi do not). 

Following LangerP^, we now identify the transition (nucleation) rate with the total 



probability current flowing across the saddle point along the unstable direction (since the 
initial density matrix has been normalized to one), i.e, the rate is given by 



oo r)n rh rh f\ r<x> 



The integral over the I ^ —1 modes (path integral) corresponds to the trace of the time 
dependent density matrix in the reduced functional space perpendicular to the unstable 
mode. 

Although we can formally find the time evolution of the density matrix by solving the Li- 
ouville equation with the Hamiltonian to quadratic order in the sphaleron basis, the resulting 
equations of motion mix the different modes and become very difficult to solve. 

It becomes clear from the above expression that, in order to compute (PD, we need to 
compute the "reduced density matrix" 

Pr{V-i,v'-i, t) = Tn^-i (e-i^V(O)e^^*) (42) 

Because in the quadratic approximation around the sphaleron solution, the Hamiltonian is 
a sum of mutually commuting Hamiltonians for each Z, it is clear that the above trace and 
thus the reduced density matrix evolves in time only through H_i 

p^it) = e-^^-^Vr(0)e^^-i* (43) 

In terms of the variable rji = (pi — (pi, the reduced density matrix initially is 

/oo _ 
^I[i^^idr]ip{7]^i,7]i;r]i^,7]i) (44) 

The time evolution of the reduced density matrix will be solely determined by the Liouville 
equation with the Hamiltonian corresponding to the unstable mode. Clearly, this is a con- 
sequence of the quadratic approximation. Thus, by taking the trace (integrating) over all 
functional directions perpendicular to the unstable mode and obtaining a reduced density 
matrix for the unstable coordinate, we have cast the multidimensional problem in terms of 
only one quantum mechanical degree of freedom corresponding to the collective coordinate 
representing the "radius of the droplet". Now the problem becomes very similar to the quan- 
tum mechanical example studied recently [^]. It is at this stage that we recognize that by 
introducing a basis which depends on the center of mass of the sphaleron (xo), we are not 
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treating translational invariance properly. For each (xo), the chosen basis corresponds to an 
orthonormal coordinate system in functional space. Changing (xq) amounts to performing a 
linear transformation on the coordinates (j)f°^ (see appendix A), thus, by integrating over all 

\x 1 

the possible values of the coordinates (f)\ one recovers translational invariance. However, 
one of the coordinates is only integrated up to the saddle point and translational invariance 
is not explicit in the above expression for the rate. To remedy this situation, it is convenient 
at this point to perform a non-linear transformation and treat ( collective coordinate. 

This is achieved by going from the (cartesian) coordinate system {0o ■, 'Pi^o] to new (curvi- 
linear) coordinates {xo ,0i^o}- There is a Jacobian associated with this transformation. To 
leading semiclassical order it is given by (see appendix A) 



J = V^sph + 0{4>) (45) 



The collective coordinate is introduced into the trace (functional integral) as|]2^, (see 
appendix A) 

Postponing the integration over the collective coordinate to the end of the calculation, 
the trace over the I ^ —1 modes can be performed easily and we find the reduced density 
matrix for the unstable coordinate to leading semiclassical order 

PriV-i,V-i) = ^exp{-^[a{r]\+7]''_,) + 2-f7]_^r]'_,]} (47) 
^ = (7^0-1,-1 - ^Q^i^-ig (48) 

I = -(K,)_^,_,-^Q^k-'Q (49) 
Qi = {K,-K2)^i,i ; 1^0,-1 (50) 
Kij, = {K^-K^y, ■ /,/VO,-l (51) 



N = JE, 



■'sph 



Det 


k' 







Det\^ 



(52) 



with K = {Ki — -ft'2) and K is the same operator but without the row and columns of matrix 
elements with the zero mode and the unstable mode in the "sphaleron" basis. Since the 
reduced density matrix evolves in time only through Hi, the continuity equation leads to 
(restoring the integral over the collective coordinate) 

J dxoj = -j dx, J_i[0_i = 0, t] (53) 
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J-i[4>-i,t] 



-ihj ( 5 



Pr(0-l,0'_l,t)U_i=0'_i 



(54) 



The current at the saddle turns out to be translational invariant (independent of Xq) 
because it only depends on 0_i (which is independent of Xo) and the integration over the 
collective coordinate will yield to a volume factor. 

The remaining task is to find the time evolution of the reduced density matrix (H^) with 
the Hamiltonian (in the Schroedinger representation) 



^ 2 



2 J,2 



(55) 



Since the initial reduced density matrix is quadratic and so is the evolution Hamiltonian, 
we propose a gaussian ansatz for the time dependent reduced density matrix (the label { — 1} 
is not explicitly written for rj to avoid cluttering of notation) 



riit) 



+ 



-[PitMt)-p*itW{t)] 

0_i-0_i(t) ; V(t)=0'_i-0_i(t) 



(56) 
(57) 



Clearly, 0_i(t) ; P(t) are the expectation value of the coordinate and canonical momen- 
tum along the unstable direction. The Liouville equation in the Schroedinger representation 
for the (—1) (unstable coordinate) now reads 



ih 



dt 



6' 



2 



(58) 



The form of the above density matrix is dictated by hermiticity. The time dependence of 
the parameters A/'(t), 0_i(t), P(t), is obtained by comparing the left and right hand side of 
(E3) and matching the coefficients of the powers of rj ; rj'. We find the following equations 



.Af_ 

AT 

—ia 

P 

0-1 



i(a-a*) 

7(0; — a*) 

; P = P* 

P 



(59) 
(60) 
(61) 
(62) 
(63) 
(64) 
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The last two equations are the classical equations of motion for the unstable coordinate in 
an inverted harmonic oscillator. The boundary conditions are: 



P(0) = ; 0_i(O) = 0_i Af{Q) = N ; a(0) = a*(0) = a ; 7(0) = 7 



In terms of the real and imaginary parts of a we find the solution to the above equations 
to be 

Qsinh[21]t] 



ai{t) 
aR{t) 

Af{t) 
tan 5 



cosh[2Qt] - cos(25) 

1 - cos(2(y) 



a(0) 
7(0) 

•AA(O) 



cosh[2ni] - cos(2(^) 

1 - cos(25) 
cosh[2Qi] -cos(25) 

1 - c()s(2(y) 



cosh[20i] - cos(2(5) 
0_i(O)cosh[Qt] ; P{t) ^ }_i{t) 



^a|(0) - 7^(0) W 
We obtain the expression for the rate per unit volume 

m 



(65) 
(66) 
(67) 

(68) 

(69) 

(70) 



L 

m 

L 

A{t) 



-n^^i{0)A{t)^E,,jM{0) exp \ -^0?.i(O)(a«(O) + 7(0))5(t) 



1 



smh[Qt\ 



1 



l + -gtanh2[nt] 



(71) 
(72) 
(73) 

(74) 



There are two competing effects that lead to the final expression for the prefactor A{t) 
in the rate. The first corresponds to the "rolling" of the expectation value of the unstable 
coordinate (0-i) down the inverted quadratic potential, the canonical momentum of this 
mode contributes to the prefactor a term proportional to sinh[i7t]. The second contribution 
has its origin in the spread of the probability distribution function (width of the gaussian) 
associated with the growth of the unstable fiuctuations contributing typical factors sinh[2f2t]. 
The width of the reduced gaussian density matrix gives the two-point correlation function of 
the fluctuation of the unstable mode. This correlation function grows exponentially because 
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of the instability. Thus we see that the spread of the probabihty distribution is the dominant 
term in the time dependent rate. This observation is in agreement with the arguments of 
classical homogeneous nucleation in that the fluctuations are the main responsibles for the 
growth of droplets and the decay of the metastable state. 

In order to find a more compact expression for the rate we need to obtain the coefficients 
a{0) ; 7(0) ; 0_i(O). In appendix B we find a simple expression for these coefficients and 
compute them in the high temperature limit and in the "thin- wall" approximation. Finally 
our expression for the rate is 

^ = Stt(^^^yA(t)e.f{-^B(t)} (76) 

with C ~ 5.1759 (see appendix B) and where we have absorbed temperature and (H) factors 
in the kernels (making them dimensionless) and defined 

^ ^ iDetipiY' 
[Det{D")_ 

D = JJ tanh 

k 

and D" the same operator in the sphaleron basis but without the rows and columns corre- 
sponding to the zero mode and the unstable mode. The rate as a function of time is depicted 
in figure (6) for an arbitrary choice of the parameters. We see that the rate vanishes at i = 0, 
reaches a maximum and eventually falls off exponentially to zero at very long times. In the 
thin wall limit, the maximum occurs at t ~ 1/W ~ 1/m. The reason why the rate vanishes 
at i = is because the initial state needs to spread out to reach the saddle and that at 
the initial time the expectation value of the canonical momentum conjugate to the unstable 
coordinate vanishes and that the density matrix has real kernels. Thus the current vanishes 
at t = 0. 

The fact that at very large times the rate must fall off to zero is a consequence of unitary 
time evolution. The total probability is conserved. As current is passing over the saddle, the 
probability density in the metastable state is depleted and thus the current must diminish. 

This is the main result of this work, a time dependent rate that incorporates the realistic 
condition of a supercooled non-equilibrium initial state. This result is certainly very different 
from the familiar approach to decay of metastable states described as steady state processes 
in thermal equilibrium. We want to emphasize the main differences with previous work on 
this problem: our expression for the rate obtains from a real time evolution of an initial 
supercooled state (ensemble). The vanishing of the rate at early times is a consequence of 
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the zero average canonical momentum and real kernels of the initial ensemble (compatible 
with local thermodynamic equilibrium), and the vanishing of the rate at late times is an 
unavoidable consequence of unitary time evolution and the depletion of probability in the 
initial state. 

Certainly, our approximation of keeping the quadratic fluctuations around the sphaleron 
configuration is not valid at long times as the amplitudes of the fluctuations become very 
large and the form of the decay rate obtained will not be accurate at long times. However, 
unitary time evolution and conservation of probability will necessarily constrain the time 
dependent rate to fall-off to zero at long times. 

This clearly illuminates the fact that the familiar results for the decay rate correspond 
to the very particular situation of a steady state and replenishing of probability in the 
metastable state. 

In a remarkable experiment on nucleation of classical fluids under shear, Min and Gold- 
burg have recently reported a time dependent nucleation rate that is strongly dependent 
on the initial conditions (shear). The nucleation rate reported by these authors has a qualita- 
tive behavior very similar to the time dependent rate ([75|), vanishing at early times, reaching 
a maximum and falling off at late times. Although our calculation is clearly quantum me- 
chanical, we conjecture that its classical limit will offer a description of classical statistical 
systems. 



4 Particle Production: 

We are now in condition to understand the production of particles of "true vacuum" as the 
metastable phase decays into the true phase via the process of thermal activation. Consider 
the operators that create and annihilate particles of momentum k of the "true vacuum", 
that is the vacuum centered at $+. Quantizing in a volume L these are given by 



y/2fiL 
1 

V2fiL 



dxe 



dxe 



ikx 



LtiHx) - $+) 



-ikx 



$ + ) + 



U(x) 



-Mx) 



(76) 
(77) 



with uj'^ the frequencies of harmonic oscillator quanta around the "true vacuum" . Now we 
can expand the field and its canonical momentum in the "sphaleron" basis by writing 



$(x) - $+ = $(x) - $_ + A$ 
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Since we expect that the maximum contribution to the rate of particle production will arise 
from the evolution of the unstable mode, we will only keep this contribution to the creation 
and annihilation operators. The average number of "true" particles as a function of time is 
then given by 

^^'^ = Trp{t) = TM^) ^^^^ 

where al{t) ; akit) are the operators in the Heisenberg picture with the Hamiltonian (|16|). 
Using, in this approximation 

0„i(t) = 0_l(O)cosh[^]^] + ^^^^sinh[fit] 
vr-i(t) = 0-1 (t) 

and that in the initial density matrix 

(r^-i)(0) = (7r_i(0)) = 



2(a(0)+7(0)) 

h 



2(«(0)-7(0)) 

and using the results obtained in the appendices we find the leading contribution of the 
unstable evolution to the total number of particles of true vacuum produced as a function 
of time (we are neglecting time independent and subleading contributions) 



k 



1 m+ 




y m 






n 


\ m ' 






Si 

+ sm\i'[nt]sJ—;^\\ (79) 

The functions Si ; 5*2 ; S3 are determined by the structure factor of the unstable mode 
f-i{x — Xo), the unstable frequency Q and the frequencies of the quanta in the true vacuum. 

In the thin wall limit the momentum integrals that determine these structure factors are 
all dominated by a peak at /c ~ of width ^ m. In this limit m~^/m~ ^ 1 ; Q/m ^ 1, in 
this limit these three functions become pure numbers of 0{1). 

Validity of the Approximations: 

Our result for the transition rate for thermal activation relies on several approximations. 
The first one corresponds to keeping the quadratic fluctuations around the droplet ( "sphaleron" ) 
configuration. This approximation provides the leading semiclassical expression for the rate 
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and is justified at early and intermediate times, since the fluctuations of the unstable coor- 
dinate will grow typically as cosh[2Qi] at early times the amplitud for the fluctuation will 
grow and presumably higher order terms (cubic and quartic) may have to be kept to fully 
understand the long time behavior. 

The thin wall approximation is justifled in the case of a small supercooling (small energy 
density difference between the metastable and the stable state) , this condition may be relaxed 
in the case of strong supercooling. In this latter case, the "droplet" conflguration will be 
indistinguishable from a localized large amplitude fluctuation. 

The initial state condition will be justified in the case of a very rapid supercooling. Since 
the maximum of the rate occurs at a time t ^ 1/m (in the thin wall limit) the initial 
supercooled state condition will be justified if the time during which the supercooling takes 
place is much smaller than 1/m. Clearly this situation is not general and will have to be 
understood case by case. In particular in the case of inflationary cosmologies, this inequality 
results in that the Hubble constant must be much larger than the mass of the quanta in the 
metastable state. This seems to be the case in the most popular theories on inflation that 
require a supercooled first order phase transition. 

Clearly at low temperatures (much smaller than Esph) tunneling will be the most impor- 
tant mechanism for metastable decay. Eventually there will be a crossover between tunneling 
and thermal activation that must be understood better. We must say however that within 
the gaussian approximation for the fluctuations around the splialcron, our solution for the 
time evolution of the density matrix does incorporate quantum corrections as may be seen 
from our result prior to taking the high temperature (classical) limit. 

Although we presented the analysis in a H- 1 dimensional fleld theory, we see no problem 
in considering our approach in 3 -|- 1 dimensions. We are currently extending our results to 
the three dimensional case and expect to report on it in a forthcoming article. 
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5 Conclusions and Implications 



In this article we have studied the process of decay of a mctastablc state via thermal activa- 
tion mediated by "sphaleron" (droplet) configurations. We offered a real time analysis based 
on the time evolution of an initially prepared density matrix corresponding to a supercooled 
state. The initial density matrix is assumed to describe the harmonic oscillations of the 
metastable state at an initial temperature much larger than the mass of the quanta in this 
state but smaller than the energy of the sphaleron configuration. 

We obtain a time dependent nucleation rate. This rate vanishes at early times as a 
consequence of the quantum mechanical spreading of the density matrix, reaches a maximum 
at a time t ^ 1/m with m the mass of the quanta in the metastable state and vanishes at 
long times as a consequence of unitary time evolution. The most important contribution to 
the rate arises from the fluctuations along the unstable direction in functional space. 

This behavior is similar to nucleation rates obtained experimentally recently in classical 
fluids under shear, but strikingly different from the familiar expressions for the decay rate. 
The difference with the usual result is that it corresponds to a steady state process in thermal 
equilibrium in which the metastable state is replenished at the same rate at which probability 
is flowing across the saddle. 

We also provided an estimate for the number of particles of the stable state produced as 
a function of time during the process of decay of the metastable state. 

For applications in the early universe or to describe processes that involve baryon number 
violation one must understand better the physics of the initial conditions. We believe that a 
steady state assumption, leading to the familiar result for the rate is not warranted in either 
case and considering an initially supercooled state may be closer to the physical situation. 

Within the context of the inflationary scenario, our results offer a rather pessimistic 
outlook. It is widely accepted that in order to complete the phase transition a fairly large 
nucleation rate is required at long times, precisely when unitarity forces the rate to vanish. 
The actual time at which the rate starts to fall off as as a consequence of unitarity will depend 
on the details of the potentials, and the amount of supercooling. A deeper understanding of 
this time regime will require to go beyond the quadratic approximation. 

With respect to sphaleron mediated baryon number violating processes, a time dependent 
rate may have important consequences for obtaining a net baryon asymmetry. In particular, 
when the thermal activation rate falls below the expansion rate in an expanding cosmology, 
sphaleron transitions go out of equilibrium. If furthermore the CP violating effects may 
be accomodated in the theory, the necessary elements for baryon number violation are in 
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place, but the question must be addressed within the realm of a realistic gauge theory with 
fermions. The influence of fermions on the dynamics of the droplets has not received much 
attention but will perhaps be an important ingredient for a deeper understanding of these 
processes. 
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6 Appendix A: 



Here we collect some relevant results for the treatment of collective coordinates. From the 
expansion around the sphaleron solution centered at Xo (0) one finds 



dxfi{x - Xo) - ^sph{x - Xo)) 

d^{x) 



(80) 



dxfi{x - Xo) 



dx 



Es,Hko + E / dxMx - a:„)^^%^0!f°l (81) 



dx 



where we have used (0). Treating collective coordinate implies |p9|, 0, 0, ^ 

^ the expansion 

$(x) = <^,ph{x -Xo) + J2 <PiM^ - ^o) (82) 
and treating {xo ; (pi^o} as coordinates. In the functional integrals (traces) the change over 



to collective coordinates is done by introducing!^ ^8 



1 = I dxo6 (0f °1) J 



J 



(83) 
(84) 
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7 Appendix B: 



In this appendix we derive an expression for the coefficients a(0) ; 7(0) that enter in the 
reduced density matrix. From the relations (^,|^) we find 

^ (a(0) - 7(0)) = [dx f dyf^.ix - x,)/_i(y - x,)e-'=(^-^) x 



' ' ^ " J 27ih 

since the Fourier transform of f-i{x) is locahzed in k-space, we can use the high temperature 
expansion of the kernels to find in the high temperature limit f3m <^ 1 

i(a(0)-7(0))^^ (86) 

We also need the sum of the coefficients; it is obtained as follows: for a fixed Xo consider 
the quantity (here rj = rj^^i) 

dyf-i{x - Xo)f-i{y - Xo)m^) - ^-){^{y) - ^-)) = {{vio?)) 

^ ^ /nz#(a:„),p(0hUh^) 

Where the expectation value is in the initial density matrix. Introducing the collective 
coordinate as in (^,^), and performing the integrals over all directions perpendicular to 
/ = — 1, this expectation value can be recast in leading semiclassical order as 

/ dx, J rf^y (r^L'r^) v.(0-i, 0-i) n 



jdxJdt^pM-iA-i) 2(«(0) + 7(0)) 

Thus 

2(a(0)+7(0)) =I'^''I'^'I ^^-'^^ - - ^°^ 2u;,tanh[/j;.l/2] ^^'^ 

Again, the Fourier transform may be calculated in the high temperature limit because the 
mode functions are localized in k-space; using tanh[f]huJk/2] ^ f3huJk/2, we finally find 

sW0) + W))«5^(5^j (90) 

With h{so) a function of only the radius of the droplet. In the thin-wall approximation, we 
find numerically 
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C = [ ~ 5.1759 

We also find in the thin-wall approximation that 

3 



0_i(O) ^ JEsph (91) 

m * 

3 Vt 

tan 5 fa 5^^=— (92) 
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Figure Captions: 

Figure 1: V{^) vs. $ for when $_ corresponds to: a metastable state (a), a degenerate 
state (b), and the true vacuum state (c). 

Figure 2: ^sph{x) vs x for a "thin- wall" sphaleron (droplet). 

Figure 3: 

Energy density as a function of the radius of the sphaleron S{x; s) in the metastable case 
for s < So {a); s — So (b) and s > So (c). 

Figure 4: 

Total energy of a droplet E{s) as a function of the "radius" s, for the cases in which 
represents the metastable state (a), degenerate state (b) and true vacuum state (c). (See 
figure 1). 

Figure 5 : 

Total energy of a droplet configuration E{s) as a function of s detailing the maximum at 
s = So- 

Figure 6: 

Decay rate per unit volume r(t)/L as a function of Wt for the arbitrary values CEgph/ ksT — 
2.0, fl/W = 0.01. (The constant coefficients in the expression for the rate were set to 1). 
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